Unveiling Order behind Complexity: 
Coexistence of Ferromagnetism and Bose-Einstein Condensation 



(N 
O 
O 
(N 



CD. Batista, G. Ortiz, and J.E. Gubernatis 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 8754-5 
(Received February 1, 2008) 

We present an algebraic framework for identifying the order parameter and the possible phases 
of quantum systems that is based on identifying the local dimension N of the quantum operators 
and using the SU(N) group representing the generators of generalized spin-particle mappings. We 
illustrate this for JV=3 by presenting for any spatial dimension the exact solution of the bilinear- 
biquadratic S=l quantum Heisenberg model at a high symmetry point. Through this solution we 
rigorously show that itinerant ferromagnetism and Bose-Einstein condensation may coexist. 
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Introduction. One of the fascinating things about mat- 
ter is the different phases it exhibits. One way or an- 
other, identifying these phases and the transitions be- 
tween them occupies the interests and activities of most 
condensed-matter physicists, physical chemists, and field 
theorists. What is typically used as a working principle 
is Landau's postulate of an order parameter (OP). While 
one generally knows what to do if the OP is known, Lan- 
dau's postulate gives no procedure for finding it. This 
problem is particularly complex if the symmetries of the 
system are hidden and "veil" the nature of the order. For 
example, it took a while to recognize that the staggered 
magnetization was the OP for anti-ferromagnetism. 

In this Letter, we present a simple algebraic frame- 
work for identifying OPs and the possible zero tempera- 
ture phases (ground states (GSs)) of quantum systems of 
spins, fermions, and bosons. Central to this framework 
is identifying the local dimension N of the quantum op- 
erators, which sets the order of the largest Lie group that 
can be used to describe those systems. The OP is con- 
structed from this hierarchical group, and the embedded 
subgroups are identified to establish a hierarchical clas- 
sification of the possible broken symmetry phases. 

To place this framework in more familiar terms, we re- 
call the Jordan- Wigner jij and Matsubara-Matsuda ||] 
spin-particle transformations that map S = | spins into 
spinless fermions and into hard-core (HC) bosons. In the 
latter case, this mapping demonstrated that ferromag- 
netism in a localized spin model was isomorphic to Bose- 
Einstein (BE) condensation in a HC boson model. We 
note that these two venerable transformations are spe- 
cial cases of the generalized transformations we recently 
introduced ||,|| and in fact they should be viewed as 
the same transformation: Their local dimension is N=2, 
corresponding to the up and down states of a spin or the 
occupied and unoccupied states of a spinless fermion or 
a HC boson. The hierarchical group is SU(2), and the 
broken symmetry states break either this symmetry in 
the ferromagnetic (FM) case or the U(l) symmetry of its 
subgroup in the BE condensation case. 

The key to the general framework we are proposing 



is the existence of a general set of SU (N) spin-particle 
transformations , fundamental for devising a general 
theoretical scheme to understand the order hidden in 
complex behavior. More simply put, through the gen- 
eralized spin-particle transformations, features that are 
subtle and hard to identify in one representation (hidden 
symmetries) can become prominent and easy to analyze 
in another (explicit symmetries). Indeed, these map- 
pings connect seemingly unrelated physical phenomena 
establishing equivalence relations among them. In the 
hierarchical group, all elements of the operator basis are 
symmetry generators || . This allows one to study the co- 
existence and competition of phases, like magnetism and 
superconductivity, with the corresponding OPs derived 
from the subgroup generators embedded in the largest 
global symmetry group of the problem. 

The N=2 case is too simple to illustrate all these 
points. We can, however, illustrate them by moving 
up to N=3 and be very concrete by studying a 5=1 
bilinear-biquadratic quantum spin Hamiltonian. In fact, 
we will present the exact solution to this problem at a 
point of high symmetry in any spatial dimension. Using 
our generalized spin-particle transformations, we show 
this model happens to map onto a lattice gas model 
of interacting particles with an internal quantum num- 
ber (two- flavored HC bosons). The exact solution trans- 
parently exhibits two simultaneously broken continuous 
symmetries associated to the formation of a polarized BE 
condensate. In this way, we rigorously show that these 
two quantum orderings can coexist . This novel result 
emerges naturally by identifying the underlying highest 
rank algebra (su(3)), providing the framework to deter- 
mine the possible broken symmetries and understand the 
role of symmetry-reducing perturbations. Using it, we 
will construct the complex phase diagram in a neighbor- 
hood of the model's high symmetry point. 

Ferromagnetism and BE Condensation. The Hamilto- 
nian we studied is [B| 
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where J < and summation is over bonds of a 

regular c?-dimensional lattice with N s sites. A spin 5=1 
operator Si is associated with lattice site i and locally 
satisfies the su(2) Lie algebra. 

The case H± — Hg=* can be conveniently written in a 
(s = i) HC boson representation {n^n^i = 6 aa >nj a ) || 

^ = ?t+^ 7! = ^+%* = ( 2 ) 

as an extended t-J like Hamiltonian, 

Hi = J E (*>lh° + H - c -) + 2J E Si • s j 

<i,j>,<r <i,j) 

+ 2jWl-^ + W) - (3) 

<iJ> V 7 

with sj = /ii/^ji/ a s = 5 operator (<x denoting Pauli 

matrices), and nj = +^jx^jx- This last form in turn 
can be rewritten as 

# 1 = 2J^P s (iJ) , (4) 

<U> 

where P s (i,j) = P s 2 (i,j) is the projector onto the sym- 
metric subspace (5=0,2) corresponding to the bond (i, j) 
which indicates that if one finds a state that is symmet- 
ric under the permutation of nearest neighbors ri and rj , 
then that state is the GS. 

For a system of Af HC bosons the GS is 

\MM,S,)) = (b^(bl i )"L\Q) , (5) 

(Af = W T + A/j. < N s ) with an energy E a = JN s z (z is 
the coordination of the lattice) and a total S z — a/ t~ a ^ , 
The operator \$ Qa is the k = component of i.e., 
b\^ a = -== J2j exp[ik • rj] 6j CT . The quasihole and quasi- 
particle excited states are 

\HW,S m )) =6 ka |* (^,5,)) quasihole, 
|*£(A/\5 Z )) = &L|*o(A/',&)) quasiparticle, U 

with the excitation energy of each being = 
Jz(| e* " — 1) where the sum runs over the vectors 
e v which connect a given site to its z nearest neighbors. 
In the |k| — > limit, Wk — * 0. 

Clearly the GS in Eq. || is a FM BE condensate with 
any partial spin polarization, and the form of the re- 
sult is independent of the spatial dimensionality of the 
lattice. We note that different values of 5 Z correspond 
to the different orientations of the magnetization AA as- 
sociated to the broken SU (2) spin rotational symmetry 
of the GS. We also note that the degeneracy of states 
with different number of particles Af indicates a broken 
U(l) charge symmetry (conservation of the number of 



particles) associated to the BE condensate ■ A signa- 
ture of BE condensation is the existence of off-diagonal 
long-range order (ODLRO) in the correlation function 

3W(y) = (bij>ja>)- When A/j and A/j. are both of order 
N s , there are two eigenvectors with eigenvalues of order 
N s and the condensate is thus a mixture. 

We can easily compute the magnetization AA and 
phase coherence of these various (non-normalized) de- 
generate GSs for a given density p = jf-. For example, 
in the fully polarized case, Af = A/j, M. = (5?) = p, 

and the ODLRO $ TT (ij) = (r; ^ rj), where 

e = 1/N S . Similarly, the two-particle correlation func- 
tion (A|A.) = $ TT (ij)^z£z^, where Aj = 5^ 
[ ]ll[ . Therefore the exact GS has two spontaneously bro- 
ken continuous symmetries (see Fig. [j]). 




FIG. 1. Coexistence of ferromagnetism and Bose-Einstein 
condensation. The correlation function $n(ij) = $ij is ex- 
pressed in the thermodynamic limit. The inset schemat- 
ically displays the SU(S) order parameter living in an 
eight-dimensional space with projections onto the SU(2) (F: 
ferromagnet) and (7(1) (BE: Bose-Einstein condensate) axes. 

The exact solution defines the features of the phase di- 
agram that our proposed framework must qualitatively 
admit. We will see below that both OPs (magnetiza- 
tion and phase), as promised, are embedded in an 5f7(3) 
OP. We remark that the phase coexistence in the boson 
representation maps back to a 5=1 FM phase coexisting 
with another spin phase. To see this consider the state 
\y o (Af,0)) for which (5f) = (5f) = (5f) = 0, which 
implies that it is a singlet state in the 5=1 representa- 
tion. We will show below that this other phase has a 
pure spin-nematic ordering. 

What can be said about the fermionic extended t-J 
Hamiltonian (i.e., Eq. ^| with HC bosonic operators re- 
placed by "constrained" fermions)? Obviously, the GS 
energies of both fermionic and bosonic Hamiltonians co- 
incide at p = 0,1. For < p < 1, one can show that 
i?o(fermion) > Eq (boson) = JN s z, indicating that the 
fermionic system phase segregates ||. 

The Framework Illustrated. In Ref. [pi, we generalized 
the spin-particle transformations to spins satisfying the 
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SU(N) Lie bracket relations. In particular, the funda- 
mental ( "quark" ) representations of SU (N) are mapped 
onto an algebra of HC bosons f| with Nf — N — 1 flavors 



Nf 



S a0 (j) = bl'b 



JV f 



N 



(7) 



where 1 < a, < iV/ run over the set of bosonic flavors, j 
is the site index, and b\ a = b\ a U (l - n i0 ). S^(j) (0 < 
v < Nf) are the components of the SU(N) spin (N 2 — 1 
linearly independent components). These components 
are generators of an su(N) Lie algebra with commutation 
relations [S^\i),S^_k)} = S^S^' (j) - S^S"*' (j). In 
particular, for N=3 



f ! - n i 1 
bt n jT 



(8) 



and up to an irrelevant constant we can rewrite Hg as 



Hg = JV2J2 [c° s e S" v (i)S vl *(X) 

<iJ> 

+ (sin e -cos e) s^{i)S^{i) 



(9) 



condensation becomes more evident: In the bosonic lan- 
guage both OPs correspond to different components of 
the 517(3) OP (see Eq. |and Fig. Q): b\ = 5 01 , b\ = S 02 , 
6| = S w , and 6| = S 20 , are the components of the OPs 
for the t and j. BE condensates, while s z — ^(S 11 —S 22 ), 
s x = i(5 12 +5 21 ), and = ±{S 12 - S 21 ), are the 
components of the FM OPs for the bosons. The ability 
of (fij) = § - (S 00 (j)) to take any value during SU(3) 
rotations of the GSs is another manifestation of the BE 
condensation. 

We will now work to answer the question: What is the 
OP of the model in the original 5*=1 language (Eq. [j])? 
We start by writing the relation between the components 
of S» v and the S=l generators of SU(2). From Eqs. | 
and H we have 

s* = JL (s 01 + s 2Q + s 02 + s w ) , 

V2 



S» 



1 



(S 



x\2 



V2i 

S z = S 11 - S 22 
2 1 



(S 01 + S w - S 02 - S w ) 



(S 12 + S 21 +S 0Q ) , (S z ) 2 -=- = 



{S x ,S y } = i{S 21 -S 12 ) , 

{S x , S z } = ±={S m - S 20 - S 02 + S 10 ) , 

{ S y, s z } = -^(s 01 - s 20 + s 02 - s 10 ) 



(ii) 



where repeated Greek superscripts are summed and 



S(j) = 



/2 - fi- 

I 3 Atj 

-hi 

\~hi 



-Hi 




(10) 



generates the conjugate representation on sublattice j 
. This expression for Hg illustrates the very impor- 
tant result that any nonlinear interaction in the original 
representation is simply a bilinear term in the new repre- 
sentation when mapped onto the highest rank algebra || . 
In particular, there are certain special points in param- 
eter space where Hg is highly symmetric. For example, 
for 6 = ? and Hg is explicitly invariant under uni- 
form SU(3) transformations on the spins [Tq] , while for 
6 = ±|£, Hg is explicitly invariant under staggered con- 
jugate rotations of the two sublattices Jig] , These sym- 
metries are hard to identify in the original spin represen- 
tation but are evident in the HC bosonic representation. 

The case where 9 = ? corresponds to the FM SU(3) 
Heisenberg Hamiltonian, and therefore the GS is the 
state with maximum total SU(3) spin S. The OP as- 
sociated with this broken symmetry is the total SU(3) 
magnetization S^ v — ^ 5 M</ (j) which has eight in- 
dependent components. When (S^) 7^ 0, the system 
orders, and the coexistence of a FM phase and a BE 



The first three operators correspond to the S=l FM OP, 
while the second five are the components of the spin- 
nematic OP (components of the bilinear symmetric trace- 
less tensor). The traceless condition of implies that 
(Sy) 2 - I = -((S z ) 2 - |) - ((S x ) 2 - |). This classi- 
fication is independent of the parity and time-reversal 
properties of A ab {t) = S a (t)S b (0) (a,b = x,y,z) 0. We 
see that by rotating the SU(3) FM OP, we find a FM or 
a spin-nematic GS for the S=l Hamiltonian. Moreover, 
we also see that the existence of a spin-nematic quantum 
phase is directly related to the BE condensation of the 
two-flavored HC bosons. Concomitantly, we have rigor- 
ously established the existence of a spin-nematic phase, 
an issue that had been outstanding in the literature jl8) . 

Unveiling Order and Establishing Complexity. With 
the phases of Hg established at the point 8 — j we will 
now detail the nature of the phases near it as revealed by 
adding a symmetry breaking field to the model. A vast 
array of rich behavior is discovered, and we now highlight 
some of the more interesting results. 

Coupling our two-flavored bosonic system Hi to a mag- 
netic field B through a Zeeman term removes the de- 
generancy of the FM GS associated to the SU(2) sym- 
metry. The resulting GS is a BE condensate polarized 
along the magnetic field direction (chosen to be the lon- 
gitudinal z direction). In fact, in the limit |B| — > 00, the 
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resulting Hamiltonian is H t -v = S(ij)(^it^jT + ^jf fyf ) 
^S(ij)("-iT — |)("jT — 5) anc ^ represents a gas of HC 
bosons with hopping t — J and nearest-neighbor density- 
density interaction V — 2 J. 

By means of the Matsubara-Matsuda transformation 
we can map H t -v onto the FM S = \ Heisenberg 
model which has an exact solution. Therefore, B removes 
the SU(3) invariance of Hi and leaves an SU(2) invari- 
ance associated to the charge degrees of freedom. In this 
way we recover the well-known connection between a FM 
S = \ state in the transverse direction and a BE conden- 
sate of spinless HC bosons. Because of the SU (2) invari- 
ance of the Heisenberg model, the FM states polarized 
along z are also GSs, so if we infinitesimally increase the 
attractive interaction between bosons, the model phase 
segregates (Ising-like ferromagnctism). 

For V = 0, H t _v reduces to an XY model whose GS is 
again a BE condensate (i.e., has off-diagonal long range 
order jp]] ) . Regardless of the dimension of the system, 
the GS of H t -v is a BE condensate for -2\t\ < V < 
and phase segregates when V < — 2\t\. In fact, for a very 
small repulsive interaction the GS is still a BE condensate 
fpif . In a similar fashion Hi is at the point separating the 
BE condensate from phase segregation, i.e., if we make 
the two-body nearest-neighbor interaction infinitesimally 
more attractive, there is phase segregation. 

A possible realization of the single-flavor bosonic sys- 
tem H t -v is the attractive U Hubbard model in the large 
\U\ limit |0| with t = -2t 2 /\U\ and V x = At 2 /\U\ (t is 
the hopping integral of the original Hubbard model). It 
is clear that if we add a nearest-neighbor attractive inter- 
action v = —2t 2 /\U\ to the original Hubbard model, we 
can exactly compute the superconducting GS and quasi- 
particle excitations in any spatial dimension and for any 
concentration of particles since V = Vi+4v = 2t. In the 
case where v = —t 2 /\U\ and V = 0, the Hamiltonian is 
an XY model whose GS displays long-range order |p| . 

Conclusions. We presented an algebraic framework 
aimed at uncovering the order behind the potential multi- 
plicity of complex phases in interacting quantum systems 
fHf . In this framework the local Hilbert space of dimen- 
sion N admits a convenient hierarchical representation in 
terms of the generators of an SU(N) group §. We illus- 
trated this for 7V=3 by first presenting the exact solution 
of a non-trivial model of interacting quantum spins. A 
key point is that from this hierarchical group (SU(3) in 
the present case), one can identify the embedded sub- 
groups (an invariant SU(2) for spin and an SU(2) for 
charge and two others connecting them), thereby estab- 
lishing a hierarchical classification of order parameters. 

Through the spin model, we rigorously demonstrated 
the possible coexistence of ferromagnetism and BE con- 
densation of s = \ hard-core bosons. This condensate 
was shown to be related to a quantum spin-nematic or- 
dering in the original spin model. These results may shed 



light on the physics associated with the recent surprising 
experimental observation of the coexistence of supercon- 
ductivity and ferromagnetism in UGe2 p3]. 

We close with the observation that a fundamental 
relation between magnetism and BE condensation has 
emerged as a consequence of our spin-particle transfor- 
mations: If a solution of a particular localized quantum 
spin S model displays long-range order in the transverse 
direction (lim| r ._ r .|^ 00 (S' i a Sj I ) ^ 0, with a = x,y), then 
the corresponding itinerant bosonic model displays long- 
range phase coherence (i.e., a BE condensation). This 
result is general and independent of the sign of the ex- 
change couplings in the spin model. 
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